Electron-magnon coupling and nonlinear tunneling transport in magnetic 

nanoparticles 
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We present a theory of single-electron tunneling transport through a ferromagnetic nanoparticle in 
which particle-hole excitations are coupled to spin collective modes. The model employed to describe 
the interaction between quasiparticles and collective excitations captures the salient features of a 
recent microscopic study. Our analysis of nonlinear quantum transport in the regime of weak 
coupling to the external electrodes is based on a rate-equation formalism for the nonequilibrium 
occupation probability of the nanoparticle many-body states. For strong electron-boson coupling, 
we find that the tunneling conductance as a function of bias voltage is characterized by a large and 
dense set of resonances. Their magnetic field dependence in the large-field regime is linear, with 
slopes of the same sign. Both features are in agreement with recent tunneling experiments. 



Metallic nanoparticles are among the best physical re- 
alizations of the concept of Fermi liquid introduced by 
Landau more than fifty years ago. Their discrete low- 
energy spectra can be put in a one-to-one correspon- 
dence with those of corresponding noninteracting elec- 
tron systems. Single-electron tunneling spectroscopy [lj 
in normal-metal nanograins provides a vivid example of 
Landau's enormous simplification of interacting Fermi 
systems. Most of the interesting phenomena studied in 
these experiments can indeed be understood in terms 
of the quantum mechanics of confined noninteracting 
quasiparticles. If the grain is made of a ferromagnetic 
transition-metal material, however, the discrete resonant 
spectrum seen in tunneling experiments 0, H| is far more 
complex than the one predicted in an independent parti- 
cle picture, and indicates that the quasiparticle states are 
coupled to the collective magnetic moment of the grain. 
Since ferromagnetic transition metals, in addition to Lan- 
dau's particle-hole (p-h) excitations, support low-energy 
collective spin excitations, it is reasonable to assume that 
tunneling transport through ferromagnetic nanoparticles 
involves some kind of spin excitations that are the finite- 
system analogue of the familiar spin-waves or magnons 
of bulk ferromagnets. So far attempts of including spin 
collective modes in tunneling transport based on a sim- 
ple toy model[l,0,Q have explained only in part the rich 
phenomena seen in experiment. 

In this paper we present a theoretical study of single- 
electron tunneling transport through a ferromagnetic 
metal nanoparticle based on a model that captures the 
salient features of its elementary excitations - p-h and 
spin collective - as derived from a recent microscopic 
study 0]. A few remarkable features seen in experiment 
emerge in a very transparent and direct way from our 
treatment of the electron-magnon coupling. We find that 
when a low-energy p-h excitation is strongly coupled to 
one of the spin collective modes, the tunneling differ- 
ential conductance versus bias voltage displays an en- 
hanced density of resonances with spacings smaller than 



the independent-electron energy mean-level spacing 5. 
The dependence of the tunneling resonances on external 
magnetic field is regulated by the behavior of the underly- 
ing quasiparticle states; it is characterized by mesoscopic 
fluctuations at small fields and a monotonic dependence 
at fields larger than the switching field. The model fur- 
ther predicts that in the limit of ultrasmall nanoparti- 
cles, where 6 is much larger than the typical magnon 
energy, the conductance should display clusters of reso- 
nances separated by an energy of order 8. 

The choice of our model is motivated by the micro- 
scopic analysis of Ref. p}, where the explicit derivation 
of the exchange-field-fluctuation propagator allows one 
to determine the elementary spin excitations (Stoner p-h 
and collective) of a magnetic grain. One finds that for 
a small nanoparticle there is one isolated spin collective 
mode below the lowest p-h excitation energy, which cor- 
responds to the ferromagnetic resonance excitation (spa- 
tially uniform q = spin- wave), of energy E ICS ~ mag- 
netic anisotropy energy/atom w 0.1 meV in cobalt. For 
large nanoparticles, the ferromagnetic resonance lies in a 
region of p-h quasicontinuum and acquires a line-width 
aE Tes , where a < 1 is the Gilbert damping parameter. 
The crossover between these two regimes occurs when 
one p-h excitation contributes to the resonance, namely 
when S = y/aE re& . Although the nanoparticles investi- 
gated in Ref. 0, Q are too small to strictly satisfy this 
condition, interactions between p-h excitations and spin 
wave modes, including the nonuniform ones {q ^ 0), will 
frequently occur. As long as the mechanism of the inter- 
action of one p-h excitation with one spin-wave mode is 
independent on the latter being uniform or nonuniform, 
we can illustrate it by following Ref. [2j|, where the uni- 
form case was considered. It was shown that when only 
one p-h excitation of energy e a b = — £ a is close to E les , 
the exchange-field propagator has two poles at energies 

El21±±± ± [ [{Elcs - e ab )/2] 
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The avoided crossing gap 7 resulting from the collec- 
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tive mode p-h coupling is found to be 7 ~ aE rcs = 
^\{b\S x \a)\, where 2S = N } - N t is the total spin 
of the nanoparticle, and Amf is the amplitude of the 
spin-splitting field. The presence of the matrix element 
I (b\ S x I o,) I in the expression for 7 emphasizes the fact that 
the coupling between spin-waves or magnon and electrons 
is ultimately due to the exchange interaction, which con- 
serves spin. Thus the quasiparticle states a and b should 
have opposite spins. In spite of the fact that most of the 
states lying close to the Fermi level have minority-spin 
character, because of spin-orbit coupling, the quasiparti- 
cle states I J) are in fact linear combinations of spin-up 
and spin-down components |7) = aj| |) |), and the 
matrix element |(6|S , 2; |a)| will not vanish. Notice, how- 
ever, that since spin-orbit interaction is relatively weak 
in transition-metal ferromagnets, the quasiparticle states 
can still be assumed to have in general one predominant 
spin character with just a small admixture of the other. 

The Hamiltonian describing the isolated nanoparticle 
in which a magnon is coupled to one p-h excitation is 

+ 1 (clc b ^ + clc a P) + Un(h-l), 

(2) 

c\ and c i with i — a,b are Fermi operators cre- 
ating and annihilating two electronic levels of energy e a 
and eb respectively, with e a < £(,. The Bose operators 
/?' and (3 describe a magnon of energy u>. Below we 
measure all energies in units of the mean-level spacing 
5 = £b — e a . The term 7 (c^c^ + c[c a j3) represents the 
electron-magnon coupling. It can be interpreted as a ver- 
tex describing an electron scattering from the electronic 
state a (respectively, b) to the state b (a) , while absorbing 
(emitting) a magnon. We will view the coupling strength 
7 as a phenomenological parameter; 7^0; represents 
strong coupling. Recently electron-boson interactions 
have been used extensively to model elcctron-phonon 
coupling in molecular single-electron transistors [8[. An 
interaction term more similar to ours has been used 
in studying magnon-assisted transport in ferromagnetic 
tunneling junctions[9J]. The last term in Eq. [^represents 
a Coulomb repulsion energy, which is nonzero when both 
electronic levels are occupied, (n) — 2. The model in 
Eq. [2l representing a double-level system coupled to one 
boson mode, is well known in quantum optics and cav- 
ity quantum electrodynamics under the name of Jaynes- 
Cummings model 10]. The model can be solved ex- 
actly, since it conserves both the number of electrons 
n = n a + nf, and the quantity (rib — n a )/2 + m, where m 
is the number of bosons. In the trivial cases n = and 
n = 2 the energy spectrum is ej^ = Lum+n/2(e a +eb+^U); 
the corresponding eigenstates are \Q,m) = (/3^) m |0) and 
|2,m) = cjjcj(/3^) m |0), where |0) is the vacuum. The 
n = 1, (rib — n a )/2 + rn + 1/2 = fc + 1 eigenspace is 
spanned by the states: 



subspace, yielding the eigenvalues e k and e k 



|l o ,A + l>=4(0t)*+i|o), \l b ,k) = cl(tf) k \0). (3) 
The Hamiltonian is now diagonalized within each k- 



4 = e° k + e Vf ±-y/ t * a + *y*{k + l), (4) 

where e res = (e b —e a )—to and e av = |(e„ + e b + lu). The 
corresponding eigenvectors are 

|±,fc) = ^ 1 ± (fc)|l a ,fc+l)+ ( 5 2 ± (fc)|l b ,fc) , (5) 



where 
5f(k) 

6t(k) = 



7\/fe + 1 



0e±-e a -^(fc + l)P+72(fc + l) 

[4-e a -u(k + l)\ 
c ±- Co -a;(fc + l)] 2 +7 2 (fc + l) 



(6) 



(7) 



On top of these states |±, k) there is also the state |1 Q , 0) 
with energy e Q , which forms a decoupled one-dimensional 
subspace in the n = 1 sector. We now assume that 
the magnetic grain is weakly coupled to metallic exter- 
nal electrodes and investigate single-electron tunneling 
transport through the grain [l[. The total Hamiltonian 
describing the system is H = Hd + Hi + H r + H t , where 
Hd is given in Eq. O Hi and H r describe the left and 
right electrodes, assumed to be normal Fermi liquids 



pa ^pa t-pce i 



a = l,r, where p is the quantum 
number specifying a quasiparticle of energy £ pQ mea- 
sured with respect to the chemical potential of lead a; 
H t is the tunneling Hamiltonian coupling the grain to 



pa^pa 



(c a + c b ) +h.c. 



Ill 



the electrodes H t = J2 P , a =r,i 
the limit of weak coupling, transport takes place via se- 
quential tunneling, which can be described by means of 
a standard rate-equation formalism for the occupation 
probabilities of the grain many-body states [ll[ . We are 
interested in the regime where Coulomb blockade is first 
lifted by applying an external bias voltage, and only the 
two charge states n = 0, 1 are involved. The master 
equations describing the kinetics of the nonequilibrium 
occupation probabilities Pj} — {P° , Pq , P^ , P^} for the 
states {|0,fc),|l«,0),|-,fc),|+,fc), fc = 0, !,•••}, 



are 



n 



pa 




E [ 2P °( W o%+k' + W %_ k , + W %) 

k' \a 

-P£W? k ,. ok - PyW« k ,. M - PSWZ fl ] ,(8) 

^'E^-^E^ (9) 



(10) 



The coefficients W a appearing in Eqs.[5]fTUlare transition 
rates between two many-body states of the grain caused 
by electron tunneling from and to the leads. For instance, 
WQ k . +k , is the transition rate from state |0, k) to \+,k') 
due to an electron tunneling from the a-electrode onto 
the grain. The W a are given by Fermi's golden rule 
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W Ok-±k' - r 0fe;±fc' n F^k' ~ e fe ~ Ma) 

x[df(k') 2 S k ,, k ^+5i(k') 2 5y, k ] , (11) 
^0% = r« a n F (e Q - MQ ), (12) 

where [i a is the electrochemical potential of lead a, which 
we assume to be shifted symmetrically around zero by the 
applied bias voltage V: fi\ = — /i r = V/2. The transition 
rates W± k ,. ok and W". Q are obtained from W^. ±fc , and 
Wg. a respectively, by replacing the Fermi function with 
[1— np] evaluated at the same energy. The tunneling rates 
r ofc;±fe' = TT E P Ka\ 2 S(e pa - (e± - e fe )) will be taken for 
simplicity to be independent of energy and lead index, 
^o k -± k i = r. The nonequilibrium steady-state probabil- 
ity PJ/ is the solution of the matrix equation MP = 0, 
where the matrix M includes all the transition rates, and 
P is a vector of all the P k 's. The dc current through the 
left or right junction is then written as 

' = (+/-) e {E^°E(<-. + < + .) 

k k' 

r k 2.^, -k;0k' r k +k;0k' 
k' k' 

+ 2P°W^ - P a W l J;} . (13) 

The transition sequence \n = 0) — > | rz = 1) — > \n = 0) 
allows the tunneling electron to probe the coupled p-h 
spin-wave excitations of the grain, which appear as reso- 
nances in the differential conductance dl/dV as a func- 
tion of the bias voltage V. We discuss first the case where 
the p-h excitation is coupled with the uniform (q = 0) 
spin-wave mode. For the nanoparticles considered in 
Refs. 0, H| ,6^1 meV, while the energy of the uniform 
spin-wave is approximately equal to the anisotropy en- 
ergy/atom ~ 0.1 meV. In Fig. [T] we plot / and dl/dV vs 
V for the case to = 0.15, which pertains to this situation. 
The calculations are done at temperature T — 0.005(5, 
corresponding to the experimental T ss 50mK. When 
7 = u) [ FigQJa) ], three sets of peaks in the conduc- 
tance are visible. The first isolated peak occurs when 
the current starts to flow, and corresponds to the suc- 
cessive transitions |0) — > |l a ,0) — > |0) which are possible 
when fn = V/2 = e a [IB|- On further increasing V, the 
current remains constant until the next lowest charging 
state €q becomes available (at eV — 6.8<5 for this case). 
For yet larger V higher states |0, k) and |— , k) acquire 
a finite nonequilibrium occupation probability, and new 
transport channels open up. In principle, each allowed 
transition |0, k) — > |— ,k) gives a resonance at e k — e k , 
as shown in the inset of Fig. [Ua), calculated at very low 
temperature, T = 0.0015. But at T = 0.005(5 only their 
envelope is visible in the form of a small bump in the con- 
ductance centered at eV = 6.9(5. The third large peak, 
appearing at eV = 7.06 is also the envelope of many 
closely spaced resonances, caused primarily by the tran- 
sitions through the second group of charged states, |+, k), 
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FIG. 1: Current and differential conductance versus bias 
voltage V for ui = 0.18. (a) The electron- magnon coupling 
strength is 7 = 0.15; (b) 7 = 0.25. The temperature T is set 
equal 0.005 5. excent in the inset, where it is eoual to 0.001 S. 
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FIG. 2: Current and differential conductance versus bias volt- 
age V, for u> — 5. (a) 7 = 0.35; (b) 7 = 0.8<5 The temperature 
T in both cases is set equal 0.005 5. 

which become available at that energy. Although values 
of 7 > ui are not very realistic, it is instructive to study 
the limit behavior of the tunneling conductance for large 
values of the magnon-electron coupling. In Fig. |Tfb) we 
plot / and dl/dV vs. V for 7 = 2uj. We can see that a 
large 7 causes the sets of resonances of Fig. QJa) to merge 
into one cluster, whose individual peaks now start to be- 
come visible also at T = 0.005(5. Notice however, that 
the mean-level spacing between the peaks is w 0.05(5, 
in fact much smaller than the experimentally observed 
resonance spacing, 0.2 5. This leads us to conclude that 
such a large density of resonances, caused by an unreal- 
istically strong coupling to the uniform spin-wave mode, 
is not the one observed experimentally. 

We now turn to the case where the p-h excitation is 
coupled to a nonuniform spin-wave mode. The exchange 
energy of the first nonuniform mode is to ~ A(a/i?) 2 , 
where A is proportional to the exchange constant, a is 
the lattice constant and R is the nanoparticle diameter. 
For a 4-nm Co nanoparticle we find uj w 1 meV, which 
is approximately equal to 5 2}. In Fig. [2ja) we plot 
the IV characteristics for the resonant case, u> = 5, and 
two different values of 7. At small 7 we have again two 
separate sets of resonances, which are now perfectly re- 
solvable even at the experimental temperature. When 7 
is increased up to 0.8, the two sets of resonances merge 
into one cluster, as shown in Fig. E^b). The number of 
resonances in the cluster is of the order of 15, with level 
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spacing w 0.35 = 0.3(eb — e a ) = 0.3w. Such a dense 
set of resonances with spacing « 0.2 — 0.55 is one of the 
characteristic features observed experimentally in tunnel- 
ing spectroscopy of magnetic nanograins. The results of 
Fig. [2] do not depend on uj being exactly equal to 5 but 
remain valid for 5 > co, although the larger u> the larger 
is 7 that takes to go from Fig. [H^a) to Fig. HJb). For 
nanoparticles much smaller than the ones considered in 
Refs.[2|, |3j, when 8 » ui, our model predicts that the 
conductance spectrum should eventually exhibit sets of 
resonances separated by an energy sa 5 cx 1/B? . 

We finally discuss the magnetic field dependence of the 
resonance spectrum. A crucial feature of our analysis is 
based on the assumption that the two bare electronic 
states |l a ,0) and |1&,0) have predominantly minority- 
spin character. The fact that minority electrons domi- 
nate the tunneling transitions had been originally pre- 
dicted in Refs. [1, Q and was later confirmed by exper- 
iments in gated devices Q. We consider first the regime 
of small external fields, where the magnetic grain is close 
to a reversal of the magnetic moment. The electronic 
states are coupled to the moment itself, and as this moves 
under the effect of the field, the energies of the states 
will be subject to random fluctuations (l2l [l3l . Ilij . Also 
the frequency of the ferromagnetic mode can fluctuate 
strongly Q. Within our model these fluctuations will re- 
sult in a quasirandom dependence of conductance res- 
onances as a function of the field. At larger fields, af- 
ter the reversal has taken place, the situation is differ- 
ent. The grain magnetic moment will point along the 
field and the energies of the minority states |l a ,0) and 
| lb, 0) will increase linearly with the field strength B, 



with a slope given by their effective g a /t, factors, which 
are w 2 since spin-orbit coupling is weak. Similarly 
the spin-wave energy dependence can be parameterized 
by u(B) = w(0) + 0/JMB-B0]- We obtain e± - e° = 
const + \[g a + g b + g ± T(B)]n B B and ef - e° k+1 = 
cor\st+h[g a +gb— gp±T(B)]fj,BB for the resonance excita- 
tion energies, where T — yj (g a — gb — gp) 2 + const/B 2 . 
If we take <?„/(, w 2 and gp < 20, we find that the ex- 
citation energies are increasing functions of B. Thus we 
conclude that the conductance spectrum exhibits essen- 
tially a monotonic linear dependence on the field, and 
the slopes of the resonance energies have the same sign. 

In conclusion, we have proposed a model that describes 
coupled electron-magnon excitations in a ferromagnetic 
metal nanoparticle. The conductance spectrum of single- 
electron tunneling exhibits a broad and dense set of res- 
onances when the coupling is of the order of the magnon 
energy. The resonant peaks show Zeeman-shifts of the 
same sign as a function of the external field. Both fea- 
tures of the model are in agreement with experiment. We 
expect that the resonances originate from the coupling to 
nonuniform spin waves; furthermore, the tunneling spec- 
trum should break into individual clusters for ultrasmall 
particles. 
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